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ABSTRACT 


An  analytical  model 


fo 


r predicting  the  pitch  and  heave 


response  of  air-cushion-supported  vehicles  (ACV)  during  overland 
operation  is  developed.  The  effects  of  air  compressibility  and  skirt 
contact  with  the  supporting  surface  are  Included  in  the  model.  The 


resulting  analytical  model  is  used  to  predict  the  response  of  a JEFF(B) 


experimental  model  to  periodically  varying  terrain.  These  predictions 


show  excellent  correlation  with  experimental  model  data,  thereby  serving 
as  verification  of  the  accuracy  of  the  analyt i ca I model  for  a 

complete  range  of  operating  conditions  for  the  JEFF(B). 


ADMINISTRATIVE  INFORMATION 

This  Study  was  sponsored  by  the  Naval  Sea  Systems  Command  under 
Task  Area  Sl4l7,  Task  li*17^;  and  administered  by  the  Amphibious  Assault 
Landing  Craft  Program  Office,  Systems  Deveiopment  Department,  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  under  Work  Unit  Number 
1- 1 180-001*. 


INTRODUCTION 


The  development  of  analytical  models  for  the  prediction  of 
the  motion  of  air  cushion  vehicles  over  random  surfaces  has 
proved  to  be  an  extremely  difficult  problem.  Only  a few  analytical 
models  have  shown  promise  of  success  (See  References  1,2, 3. 3)  The 
difficulties  in  developing  a suitable  analytical  model  are  attributed  to 
the  strong  nonlinear  characteristics  of  the  physical  phenomena  in 
question.  These  nonlinearities  arise  from  the  deformation  of  the 
vehicle  skirts,  the  complicated  air  flow  patterns  In  the  supporting 
cushions,  and  in  general  the  nonlinear  characteristics  of  the  Interface 
between  the  craft  and  supporting  surface.  If  the  supporting  surface  is 
deformable,  additional  nonlinearities  arise.  In  order  to  eventually 
produce  craft  having  optimum  design  characteristics,  it  is  necessary  to 
develop  analytical  tools  that  accurately  model  craft  behavior.  This 
report  presents  the  successful  development  of  a nonlinear  time-domain 
model  that-  can  be  used  for  predicting  craft  response  in  pitch  and  heave 
during  overland  operations.  The  analytical  model  developed  does  not 
contain  any  artificial  damping  or  restoring  coefficients,  artificial 
inertial  or  mass  properties,  or  artificial  tuning.  All  parameters 
employed  are  derived  directly  from  the  physical  system  under  study.  The 
degree  of  success  of  the  current  analytical  model  Is  indicated  by  the 
excellent  correlation  results  achieved  with  experimental  model  data.  The 
simulations  conducted  for  purposes  of  verification  Included  prediction 
of  pitch  and  heave  for  motion  over  sinusoidal  terrain,  trapezoidal 
terrain,  ramps,  triangular  bumps,  and  step  function  changes  in  terrain. 
The  more  interesting  sinusoidal  results  are  presented  in  this  report. 


*References  are  listed  on  page  28. 
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The  analytical  model  presented  In  this  report  Is  appro- 
priate for  a model  scale  craft.  Once  the  analytical  model  Is  fully 
developed,  generalization  to  full  scale  requires  only  a straightforward 
process  of  changing  geometrical  parameters  to  full  scale  values.  The 
full  scale  analytical  model  can  be  used  in  an  extremely  cost  effective 
manner  for  design  optimization  studies,  trade-off  studies,  and 
performance  predictions. 

DEVELOPMENT  OF  THE  MATHEMATICAL  MODEL 
The  air  cushion  vehicle  configuration  considered  In  the  analytic 
model  consists  of  a fully-skirted  perimeter  with  two  supporting  air 
cushions  separated  by  a transverse  stability  skirt.  The  air  cushions 
and  peripheral  skirts  are  supplied  with  air  from  a fan  chamber,  which 
in  turn  is  supplied  with  air  from  the  lift  fan.  This  general  configur- 
ation, which  is  illustrated  in  Figure  1,  is  representative  of  the 
Amphibious  Assault  Landing  Craft  (AALC)  Program  JEFF(B)  craft. 

All  mathematical  models  are  developed  in  a body  centered 
coordinate  system  having  the  z axis  directed  upward  and  the  x axis 
forward.  Positive  pitch  corresponds  to  a bow  up  orientation.  The 
analytic  model  has  been  developed  to  compute  ACV  motion  in  three  degrees 
of  freedom,  specifically  pitch,  heave,  and  surge.  Only  pitch  and  heave 
results  are  presented  in  this  report  since  experimental  data  including 
surge  was  not  available  for  validation  purposes.  Means  have  been  developed 
to  extend  the  present  mathematical  model  to  include  all  six  degrees  of 
freedom. 
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The  basic  equation  of  motion  for  heave  Is  given  by 


mz  - + Fy  - mg  (l) 

where*  the  products  P^S^  and  represent  the  vertical  force  arising 

from  the  forward  and  aft  cushions,  represents  the  vertical  force 
arising  from  skirt  contact  with  the  supporting  surface,  and  mg  repre- 
sents the  gravitational  force.  The  differential  equation  describing 
pitch  is  given  by 

l§-  P,S,(x,-Xg)  . (Xj  - Xg)  . « 

t P,B(h^,  - h„3)  {Zg  . 4 (h^,  t h^j))  (2) 

* " ^W2^  ^^9  * ^ ^’’k3  ^k2^^ 


The  first  two  expressions  appearing  on  the  right  hand  side  of  the 
above  equation  represent  the  pitching  moment  arising  from  the 
forward  and  aft  cushion  pressure  forces  respectively.  M represents 
the  moment  produced  by  skirt  contact  with  the  supporting  surface 
and  arising  from  air  escaping  under  the  skirts.  The  last  two  terms 
represent  moments  identified  with  the  pressure  differentials  that 
exist  across  the  skirts. 


* A complete  description  of  ail  symbols  appears  in  the  notation  section. 

% 


The  above  second  order  differential  equations  are  solved 


directly  by  a Taylor  series  expansion  using  evenly  spaced  time  steps, 

At.  Extensive  sensitivity  studies  were  conducted  to  determine  an 
optimum  integration  time  step  and  to  assure  that  Integration  Instabili- 
ties were  not  encountered.  Various  exciting  forces  were  used  during 
these  sensitivity  studies.  In  each  case,  the  step  size  At  was  decreased 
until  variations  In  final  results  remained  less  than  1 percent.  A time 
step  size  of  0.0025  seconds  proved  sufficient  for  1 percent  accuracy  for 
all  cases  examined  in  a model  scale  simulation. 

Expressions  for  the  various  forces,  moments,  and  pressures 
appearing  on  the  right  hand  side  of  the  above  differential  equations  are 
developed  in  subsequent  sections.  The  expressions  are  implemented  in  a 
digital  computer  program  which  is  used  to  generate  the  predicted 
results  presented  In  this  report. 

CUSHION  PRESSURE  MODEL 

The  pressure  in  each  chamber  of  the  segmented  cushion  of  the  air 
cushion  vehicle  is  assumed  to  be  spatially  uniform,  with  forces  acting 
at  fixed  centers  of  pressure,  and  T2.  In  the  actual  ACV  there  will 
exist  a pressure  gradient  in  each  of  the  cushions  arising  mainly  from 
the  mass  transfer  of  air  (4,  5,  6).  This  pressure  gradient  will  depend  in 
part  on  craft  orientation  and  velocities  and  In  part  on  the  changing 
geometrical  character  Is t Ics  of  the  solid  surface  supporting  the  craft. 

The  changing  pressure  distribution  will  alter  the  location  of  the 
center  of  pressure  during  craft  motion. 
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The  current  approx! mat  ton  Is  adequate  provided  the  pressure 
gradients  are  small.  A detailed  examination  of  the  effects  of  the 
pressure  gradient  would  require  solving  a set  of  partial  differential 
equations  with  temporally  and  spatial  I y dependent  coefficients. 

A numerical  procedure  of  this  type  is  currently  being  investigated 

under  a separate  effort. 

For  a given  fan  air  flow  rate,  Q,  the  pressure  in  the 
intermediate  chamber  between  the  fan  and  the  forward  and  aft  cushions 
Is  assumed  to  be  represented  by  a simple  quadratic  form, 

Pf  - + CjQ  + (3) 

This  mode?  takes  into  account  ducting  effects  and  fan  characteristics. 
The  coefficients  used  in  the  current  study  were  obtained  from 

experiments  conducted  by  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC)  on  a model  of  the  JEFF(B),  an  Amphibious 
Assault  Landing  Craft  (6).  This  simple  nonlinear  equation  was  chosen 
for  convenience  and  can  be  easily  replaced  by  tabulated  values,  more 
complicated  fan  maps,  or  analytic  models  of  the  dynamic  response  of  the 
lift  fan  and  air  distribution  systems. 

The  pressures  in  the  forward  and  aft  cushions  are  obtained 
from  air  flow  continuity  equations  and  from  discharge  equations  for  the 
various  orifices  between  the  fan  chamber  and  each  of  the  cushions, 
between  each  of  the  cushions  and  the  atmosphere,  and  between  forward 
and  aft  cushions.  The  continuity  equation  specifies  that  the  sum  of  the 


input  flow  to  a cushion,  the  output  flow,  and  the  qains  or  losses  due  to 
compressibility  must  equal  zero.  The  equation  of  continuity  for  the 
forward  cushion  Is 

'J.l  - '’ll  * * '=cl  ■ “ 

and  for  the  aft  cushion 

'i.2  ■ «I2  * V - 1„2  - 5.3  ♦ «C2  ■ » 

where  Qj^  and  are  the  Inlet  flows  from  the  fan  chamber,  is  the 
stability  skirt  flow  between  the  cushions,  and  ^re  escape  flows 
through  forward  and  aft  skirts,  and  flows  associated 

with  changes  In  cushion  volume  due  to  craft  motion,  and  Q ^ are  the 
air  flows  resulting  from  volume  changes  arising  from  varying  terrain, 
and  and  Q^2  ®re  effective  air  flows  associated  with  air  compressibility 
effects  and  the  rate  of  change  of  density  of  the  cushion  fluid. 

The  relationship  between  pressure  and  air  flow  through  various 
openings  is  specified  by  the  general  orifice  equation 

0 = CpjA./2AP7p  (5) 

where  Cg.  is  an  orifice  discharge  coefficient  associated  with  orifice 
area  Aj , and  AP  is  the  pressure  difference  across  the  orifice. 

The  discharge  coefficients  used  In  the  present  analysis  have 
been  derived  from  experimental  studies ' currently  underway.  Geometrical 
data  were  obtained  directly  from  the  1:12  scale  JEFF(B)  model  used  in  the 
validation  experiments.  Orifice  areas  consisted  of  a large  number  of 
small  orifices  of  various  shapes.  The  measured  leakage  area  represented 


a lower  bound  since  small  leakage  areas  could  not  be  accurately  measured. 
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The  escape  areas  under  the  skirts  and  the  stability  seals 


are  calculated  by  Integrating  the  gap  between  the  bottom  of  the  skirt 
and  the  supporting  surface  under  the  skirt  perimeter.  The  resulting  leakage 
area  depends  upon  the  craft  motion  and  craft  orientation,  which  are 
computed  In  the  motion  simulation,  and  also  upon  the  terrain  geometry 
relative  to  the  ACV.  The  analytical  model  for  the  leakage  gap  employed 
in  the  simulation  results  presented  In  this  report  Is  given  by 

a(x,z,  t)  - z + >^)  - h|^(x)  - z^(x  + Ut)  (6) 

where  hj^  Is  the  skirt  height,  z^  Is  the  terrain  profile,  U Is  the  craft 
velocity,  and  0 is  the  pitch  angle  which  has  been  assumed  to  be  small. 

The  small  angle  approximation  Is  appropriate  since  predicted  pitch  angles 
are  typically  less  than  6 degrees.  A negative  value  for  & indicates 
that  the  skirt  is  in  contact  with  the  surface,  i.e.,  there  Is  no  gap, 
in  which  case  £ Is  set  to  zero.  In  order  to  handle  the  case  in  which 
il<0,  I can  be  written  as  where  H(£)  is  the  Heaviside  function. 

AIR  FLOW  DUE  TO  CRAFT  PUMPING 

The  craft  pitch  and  heave  motion  resulting  from  interaction  with 
Its  environment  causes  the  volume  of  the  forward  and  aft  cushions  to 
change  as  a function  of  time.  The  contribution  to  incompressible  flow 
resulting  from  craft  motion  Is  assumed  to  be  equal  to  the  product  of  the 
craft's  average  vertical  velocity  In  the  chamber  In  question  and  the  cross- 
sectional  area  of  the  forward  or  aft  cushion,  as  appropriate.  The  craft 
pumping  term  Is  therefore  given  by 

- {z  + (k  - Xg)  0}S  (7) 
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where  S is  the  cushion  cross“scctIonal  area  and  H Is  the  location  of 
the  cushion  centroid.  The  centroid  of  the  cushion  Is  assumed  to  be  at 
the  center  of  the  cushion  segment  under  examination. 

AIR  FLOW  DUE  TO  TERRAIN  PUMPING 

The  air  transport  effect  of  forward  motion  over  fixed  terrain  is 
given  for  each  cushion  segment  by 


0^  - /Z  z^(x,t)  B dx 

^ I W 


^ Li  3t 

where  B is  the  craft  beam,  x Is  measured  along  the  length  of  the  craft, 
and  Lj  and  L^  are  the  x positions  of  the  cushion  front  and  back  boundaries. 
For  free  surface  simulation  this  term  is  referred  to  as  wave  pumping  and 
must  include  the  celerity  of  the  wave  field. 

The  form  of  the  above  expression  Is  of  particular  interest 
when  the  terrain  has  a simple  trigonometric  profile  since  arbitrary 
profiles  can  always  be  represented  by  a suitable  Fourier  series.  If 
the  terrain,  is  given  by 

z^(x,t)  = a cos(k^  + wt)  (9) 

and  the  integration  occurring  in  the  expression  for  is  carried  out, 

the  terrain  pumping  term  for  the  bow  cushion  becomes 
a Bu 

= -1 — ^(cos  k L/2  - 1)  cos  U)t  - sin  k L/2  sin  U5t} 
and  for  the  aft  cushion 


Qw2  * ~r — Hi  - cos  k^L/2)  cos  wt  - sin  k^L/2  sin  »t}  (lOb) 


AIR  FLOW  DUE  TO  COMPRESSIBILITY  OR  CUSHION  PUMPING 

The  effect  of  the  compressibility  of  the  air  In  the  cushion 
of  the  ACV  can  be  represented  to  first  order  In  time  rate  of  change  of 
density  by  an  air  flow  term  given  by 

« 

- VO/p^  („) 

where  p Is  the  cushion  air  density  and  p.  Is  the  ambient  air  density. 

An  adiabatic  Isentropic  equation  of  state  Is  employed  to  relate  the 
density  to  the  pressure  In  the  cushion.  The  state  equation  Is  given  by 

1 ♦ p/Pg  - (P/Pa)^ 

where  p^  Is  the  ambient  atmospheric  pressure  and  y Is  the  ratio  of 
specific  heats.  If  the  state  equation  Is  solved  for  p/p^  the  result 
Is 

P/P-  - -i-  (£-)  *=  _L 

Yfta  Pa  Y yp. 

where  the  last  expression  Is  obtained  by  neglecting  terms  of  the  order, 


The  effective  air  flow  due  to  compressibility,  which  Is  the  true 

cushion  pumping  term  can  be  determined  as  a function  of  pressure  by  using 

the  above  expression  for  p/p  , resulting  in 

d 

Qc  = Vp/pa  Vp/(ypa)  (12) 

SKIRT  FORCES 

The  analytical  model  for  the  force  arising  from  skirt  contact 
with  the  supporting  surface  is  highly  nonlinear  even  though  simplifying 
assumptions  are  made  in  modeling  the  skirt  geometry,  it  is  assumed  that 
the  materials  used  in  constructing  the  skirts  are  sufficiently  flexible 
that  restoring  forces  associated  with  the  material  can  be  neglected.  It 
is  also  assumed  that  frictional  forces  arising  from  contact  with  the 
supporting  surface  do  not  deform  the  skirt  in  the  body  x direction, 
implying  that  when  contact  is  made  the  skirt  compresses  in  the  vertical 
direction  only.  An  illustration  of  the  skirt  deformation  during  contact 
with  the  supporting  surface  is  presented  in  Figure  2.  The  vertical  force 
associated  with  the  supporting  surface  contact  arises  from  the  cushion 
or  skirt  bag  air  pressure  acting  over  the  surface  of  contact.  In  terms 
of  the  variables  illustrated  in  Figure  2,  the  total  vertical  force 
arising  from  the  bow  skirt  contact  is  given  by 

Fj  = P,B(d2  - 6^)  + 2PfBd^ 

where  B is  the  skirt  beam,  B(d2  ~ d^)  is  the  skirt  contact  area,  2Bd^ 

is  the  skirt  bag  contact  area,  P^  is  the  cushion  pressure,  and  P^  Is 
the  skirt  bag  pressure  which  is  assumed  to  be  the  same  as  the  fan  chamber 
pressure.  The  above  expression  can  be  written  In  terms  of  the  skirt 
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height  and  the  angle  that  the  skirt  makes  with  the  vertical.  The 
resultant  force  is  given  by 

P,  = PiBdj  tan  - dj)  + ZP^Bd^  (13) 

The  distance  d^  can  be  obtained  from  the  effective  radius  of  the  bow 
skirt  bag  and  the  bow  skirt  deformation  height  hence, 

^3  -{(i'b  - L,)  (2R^  - i'b  + 

The  aft  skirt  is  modeled  without  a bag.  Hence,  the  vertical 
force  associated  with  aft  skirt  contact  is  given  by 

Fz  = (Pf  - Pj)  tan  0^) 

where  ?2  is  the  aft  cushion  pressure,  Is  the  deformation  of  the  aft 
skirt,  and  is  the  angle  the  aft  skirt  makes  with  the  vertical.  The 
total  vertical  force  Is  then  given  by 

Py  * ®b  ' ^^f*^3  ^2®*'* a "a 

Occasionally  the  craft  hull  above  the  skirt  bag  will  impact  the 

supporting  surface.  When  this  occurs,  the  vertical  force  is  replaced 
by  a rapidly  increasing  function  of  the  deformation  height.  The  defor- 
mation of  the  skirts,  2',  can  be  written  In  terms  of  the  mathematical 
expression  given  for  in  Equation  (6) 


i^'b  “ -^(L/2,  z,  t)  H(-Ji) 

- -Jl(-L/2,  z,  t)  H(-^) 

9 

where  H(i^)  is  the  Heaviside  function. 
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SKIRT  MOMENTS 


Three  different  effects  are  considered  In  modeling  the  moments 
due  to  the  skirts.  These  effects  are  the  reaction  phenomena  arising  from 
air  escaping  from  the  cushion  under  the  skirt,  frictional  effects  due 
to  skirt  contact  with  the  supporting  surface,  and  skirt  deformation 
effects  during  skirt  contact  with  the  supporting  surface.  The  reaction 
moment  Is  given  by 

Mr  - -P,Bil(L/2,  2,  t){2_  + h.  , + H(L/2,  z,  t)}  H(^) 

® (16) 

+p2B2(-L/2,  z,  t){zg  + z,  t)}  H(£) 

The  seal  contact  moment  Is  obtained  by  multiplying  the  seal  contact 

forces  by  appropriate  moment  arms  and  summing.  The  results  are 

Mj.  = F,(L/2  + tan  ‘ ^^7) 

where  Equations  (13)  and  (1l»)  have  been  employed.  The  moment  arising 
from  frictional  forces  Is  given  by 

Mp-F,p(h^,  - »'^)  - - £',)  (18) 

where  y 's  the  coefficient  of  friction.  The  total  skirt  moment  M is 
given  by 

M - M^  + M^  + M^ 


13 


SOLUTION  OF  AIR  FLOW  EQUATIONS 


.1 


; ) 


The  continuity  equations  for  the  forward  and  aft  cushions 
given  by  Equation  (k)  can  be  written  as  functions  of  the  unknown 
pressures  and  by  using  Equations  (3),  (5).  (8),  (10),  and 

(12).  The  explicit  form  for  the  resulting  pressure  dependent  functions 
Is  not  presented  here,  since  evaluations  of  the  functions  In  the 
computer  simulation  program  proved  more  efficient. 

The  three-dimensional  generalized  Newton-Raphson  numerical 
procedure  Is  used  to  solve  the  resulting  equations  for  the  pressures 
Pp  P^,  and  P^  required  In  the  equations  of  motion.  Since  the  equations 
Involved  do  not  contain  very  complicated  expressions  In  terms  of 
pressures,  a more  sophisticated  numerical  procedure  Is  not  warranted. 
Convergence  Is  occasionally  slow  due  to  oscillatory  roots.  In  these 
cases,  convergence  Is  accelerated  by  automatically  Initiating  the  Iterative 
procedure  with  a priori  values  equal  to  averages  of  the  oscillatory  roots. 
This  only  occurs  when  the  computed  values  are  close  to  the  solution. 

The  algorithm  used  In  Implementing  the  Newton-Raphson  procedure 
has  other  features  which  Improve  convergence.  The  capability  of  de- 
creasing the  predicted  refinement  by  a given  fraction  proved  quite 
useful  In  eliminating  a frequently  encountered  problem  characteristic  of 
many  Newton-Raphson  algorithms,  namely  overestimating  parameter  refinements. 
Time  histories  of  the  cushion  pressures  generated  by  the  Newton-Raphson 
solution  procedure  were  examined  In  detail  for  a wide  variety  of  cases. 

All  values  for  the  pressures  and  the  time  rate  of  change  of  the  pressures 
were  always  consistent  with  the  corresponding  craft  orientation  and 
motion  variables. 

\k 


EXPERt MENTAL  VALIDATION  OF  ANALYTIC  MODEL 


Extensive  correlation  studies  with  experimental  data  were 
carried  out  In  order  to  verify  the  analytical  model  presented  In  this 
report.  A 1/12  scale  model  of  the  JEFF(B)  craft  was  used  by  DTNSROC 
In  an  experimental  program  In  which  the  model  was  tested  over  a wide 
range  of  terrain  types  representative  of  obstacles  that  might  be  en- 
countered In  a realistic  overland  operating  environment.  The  experiments 
Included  towing  the  model  over  solid  periodic  waveforms  of  a wide  range 
of  amplitudes,  wave  lengths,  and  encounter  frequencies.  The  results  of 
this  experimental  program  are  presented  In  Reference  7.  The  model  was 
constrained  In  surge  during  the  experiments  used  for  the  correlation 
results  presented  In  this  report  and  hence  numerical  results  are  pre- 
sented only  for  pitch  and  heave  response. 

The  numerical  Implementation  of  the  analytic  model  was  exercised 
as  presented.  No  artificial  restoring  or  damping  coefficients  were 
Introduced,  a practice  that  Is  all  too  common  In  mathematical  modeling. 
Computed  response  magnitudes  were  not  modified,  and  no  attempt  was  made 
to  adjust  mean  value  responses.  The  analytic  model  was  used  In  as 
physically  pure  a form  as  could  be  developed  and  the  result  Is  that 
exceptional  agreement  Is  obtained  between  the  experimental  data  and  the 
analytic  predictions.  The  successful  verification  which  has  been  obtained 
Is  essential  as  a demonstration  of  the  value  of  the  analytic  model  In 
the  overall  effort  of  research,  design,  and  testing  of  alr-cushlon- 
supported  vehicles. 
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f (Tient  between  analytical  model  predictions  and  experimental  data  has  not 

* been  demonstrated  previously  or  presented  in  the  literature  for  air-cushion- 


sup|.K)rted  vehicles.  Discrepancies  occur  only  in  regions  where  slopes 
are  rather  steep,  making  it  difficult  to  carry  out  accurate  experimental 
measurements.  The  predicted  trends  with  respect  to  speed  agree  very 
well  with  the  measured  trends  in  these  regions,  and  in  these  and  all  other 
regions  the  predicted  results  are  well  within  the  expected  experimental 
accuracy . 

General  trends  that  are  representative  of  the  expected  craft 
behavior  are  quite  apparent  in  the  analytical  results  presented  in  this 
report.  In  the  long-wave- length  limit,  it  is  expected  that  the  craft 
motion  will  be  parallel  to  the  wave  surface.  Consequently,  both  the  pitch 
and  heave  transfer  functions  approach  unity  in  the  long-wave- length  limit. 
As  the  wavelength  approaches  small  values,  the  effects  of  the  peaks  and 
troughs  under  the  craft  tend  to  average  out,  producing  no  net  moment 
or  heave  force.  Hence  in  the  short-wave- length  limit  the  envelope  which 
contains  the  maxima  of  the  transfer  functions  should  finally  approach 
zero. 

The  detailed  behavior  of  the  response  includes  null  points  at 
integral  fractions  of  wave-length-to-craft-length  ratio.  These  null 
points  arise  from  a resonance  type  of  phenomenon  between  the  cushions 
and  the  wave  terrain,  resulting  in  cancellation  of  contributions  to  the 
force  and  moment  associated  with  bow  and  aft  cushions  and  with  null 
conditions  for  craft  and  terrain  pumping  terns  so  the  null  points  will 
be  sharp  if  the  craft  length  is  well  defined.  In  the  present  analytical 
model  for  the  skirts,  allowance  Is  made  for  craft  length  changes  due  to 
skirt  contact  with  the  terrain.  These  changes  in  length  tend  to  reduce 


17 


the  sharpness  of  the  null  points,  an  effect  that  Is  quite  apparent  In 
both  the  predicted  results  and  the  experimental  data.  This  type  of 
response  behavior  Implies  that  the  skirt  model,  although  simple.  Is 
va 1 id. 

APPLICATIONS  OF  THE  ANALYTICAL  MODEL 
An  air  cushion  vehicle  motions  or  simulation  model  may  be 
employed  with  assurance  in  craft  dynamics  studies  only  after  it  has  been 
thoroughly  validated  through  correlation  with  experimental  data.  Having 
shown  that  the  present  ACV  mathematical  model  performs  correctly  for 
overland  operation  in  two  degrees  of  freedom  over  rigid  periodic 
waveforms,  it  is  appropriate  to  present  a brief  Investigation  of  various 
dynamic  characteristics  of  the  ACV.  Included  in  this  overview  of  dynamic 
response  characteristics  is  an  Investigation  of  linearity  in  the  pitch- 
heave  response,  an  examination  of  the  magnitude  of  terrain  pumping  terms 
in  the  mathematical  model,  a demonstration  of  the  effect  of  modifications 
. in  the  properties  of  the  fan  curve  on  vertical  plane  motion,  a brief 

examination  of  the  effect  of  compressibility  of  the  cushion  fluid,  and  an 
initial  investigation  of  the  subharmonic  character  of  ACV  overland 


response. 

t 

« 

^ RESPONSE  LINEARITY 

h 

M The  degree  of  nonlinearity  in  the  overland  response  dynamics 


/ *• 
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of  the  ACV  analytical  model  has  been  Investigated  as  a function  of  terrain 
amplitude  for  sinusoidal  excitation  over  a range  of  encounter  frequencies. 
The  results  of  this  Investigation  are  presented  in  Figures  9 through  12. 
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Pitch  and  heave  transfer  functions  are  shown  In  Figures  9 and  10  as 
functions  of  wave  slope  for  fixed  values  of  wavelength.  Results  are 
obtained  for  various  frequencies  by  changing  the  speed  of  transit  over 
the  wave  field.  For  a linear  system  the  transfer  functions  would  appear 
as  straight  horizontal  lines,  Indicating  that  the  pitch  and  heave  response 
at  a given  encounter  frequency  Is  directly  proportional  to  the  wave 
amplitude.  The  deviations  from  constant  response  appearing  In  these 
figures  Indicate  the  degree  of  nonlinearity  observed  at  specific  excitation 
frequencies.  It  Is  evident  that  for  any  specific  excitation  frequency 
there  are  regimes  In  which  the  ACV  response  may  be  classified  as  linear 
and  other  regions  In  the  wave  amplitude  domain  in  which  the  nonlinear 
behavior  of  the  vehicle  Is  dramatic.  Factor-of-two  variations  In  the 
response  amplitude  are  not  uncommon  for  small  wave  slopes,  but  the  ACV 
response  tends  to  be  more  linear  for  steeper  waves.  It  should  be  noted 
that  the  peaks  in  the  pitch  transfer  functions  shift  toward  larger  wave 
amplitudes  as  the  craft  speed  Is  decreased.  It  Is  especially  Important 
to  note  the  manner  In  which  the  heave  Inflection  point  moves  to  larger 
wave  slopes  with  decreasing  speed.  The  curve  presented  for  A6.8  knot 
operation  (13.5  kts  model  scale,  F^*1.56)  shows  that  the  magnitude  of 
the  heave  response  doubles  In  the  middle  of  a range  of  wave  slopes 

representative  of  common  terrain.  ! 

The  predicted  transfer  functions  are  presented  as  functions  of 
encounter  frequency  In  Figures  11  and  12  for  three  different  wave 

I 

amplitudes.  In  this  presentation,  all  three  normalized  curves  would  | 

coincide  for  linear  system  response.  The  large  differences  among  these  | 
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curves  indicate  stronqly  nonlinear  behavior  exhibited  by  the  ACV  as  the 
wave  amplitude  approaches  zero.  Over  the  range  of  amplitude  ratios  (a/L) 
considered,  the  frequency  associated  with  the  pitch  resonance  shifts  to 
hiqher  frequencies  for  flatter  terrain  (i.e.,  smaller  wave  amplitudes). 

The  peak  heave  resonance  frequency  behaves  in  a similar  manner,  moving 
to  a hiqher  frequency  with  decreasing  wave  slope.  The  implication  of  this 
behavior,  in  which  the  nonlinear  character  In  ACV  response  dominates 
the  dynamic  performance  for  the  flatter  terrain,  is  disturbing  since  it 
is  under  these  conditions  that  one  expects  linear  theory  to  be  applicable. 
Fortunately  the  concern  Is  only  academic  since  the  actual  ACV  response 
is  minimal . 

TERRAIN  PUMPING 

The  expression  "terrain  pumping"  refers  to  changes  in  the  Cushion 
air  flow  due  to  motion  of  the  craft  over  terrain  of  varying  elevation. 

For  a sinusoidal  terrain,  the  terrain  pumping  effect  results  in  a sinu- 
soidal air  flow  term  (Q^)  in  each  cushion.  When  the  bow  and  aft  cushion 
discharges  are  in  phase  or  additive,  the  heave  response  is  magnified. 
Likewise,  whenever  the  changes  in  bow  and  aft  cushion  are  out  of  phase, 
the  pitch  response  will  be  most  affected.  The  time-average  of  the  sum 
and  difference  of  the  two  flows  is  shown  in  Figure  13.  Examination  of 
the  flows  shows  that  the  difference  peaks  at  a wave-to-cush ion- length 
ratio  of  1 (X/L*l)  and  has  a null  point  at  the  half  craft  length  (>/L“0.5) 
while  the  sum  peaks  at  a wave-to-cushion- length  ratio  of  2 (X/L“2)  and 
has  a null  point  at  unity  (X/L«1).  This  phenomenon  explains  the  appearance 
of  null  points  in  the  pitch  and  heave  response  but  since  both  pitch  and 
heave  are  strongly  coupled,  the  maximum  responses  are  expected  to  occur  at 
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wave-to-cush!on-length  ratios  shifted  somewhat  from  these  Integer  values. 
This  Is  quite  evident  from  the  transfer  function  results  presented  In 
Figures  ^ through  8. 

Other  factors  also  determine  where  the  maxima  and  minima  will 

occur.  Among  these  are  the  resonant  response  of  the  craft  In  pitch  and 

heave  and  the  Interaction  of  the  skirts  with  the  supporting  surface.  This 

Interaction  results  In  a speed-dependent  minimum  In  the  heave  response 

at  a wave-to-cush Ion-length  ratio  of  1.3  for  a full-scale  JEFF(B)  speed 

of  27.7  kts  (F^*0.93)  and  at  about  1.1  for  a speed  of  52  kts  (F^*1.7^). 

This  result  may  be  Interpreted  as  follows.  The  contact  forces  given  by 

Equations  (13).  (1^),  (I6).  (17),  and  (I8)  are  highly  nonlinear  and  thus 

can  be  expected  to  shift  the  minimum  response  point  away  from  unity. 

Since  the  wave  pumping  term  Increases  linearly  with  speed  or  encounter 

frequency,  It  would  tend  to  dominate  the  response  at  higher  speeds. 

It  Is  observed  In  Figures  6 and  8 that  the  heave  minimum  Is  closer  to 

a reduced  frequency  of  1 at  52  kts  (F  =1.7^)  than  It  is  at  27.7  kts 

n 

(Fp-0.93). 

LIFT  FAN  CHARACTERISTICS 

It  is  Important  In  the  design  of  an  ACV  lift  system  to  determine 
the  effects  that  various  fan  characteristic  curves  will  have  on  the 
dynamic  performance  of  an  ACV.  The  present  ACV  analytical  model  may  be 
employed  to  parametrically  investigate  lift  system  and  fan  flow  modifi- 
cations. Two  numerical  experiments  have  been  performed:  changes  In  the 

total  flow  (()))  at  a given  pressure  (i|))  and  changes  In  the  slope  (d<li/d<l>) 
of  the  pressure-flow  curve  at  the  mean  operating  point.  The  fan  maps  used 
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Jn  these  simulations  are  shown  in  Figure  1l*  and  the  dynamic  response 
characteristics  associated  with  the  fan  variations  are  shown  In  Figures 
15  and  16.  The  response  characteristics  should  be  compared  with  the 
response  resulting  from  the  nominal  fan  performance  curve.  The  major 
change  in  the  pitch  response  occurs  for  a decrease  in  the  lift  system 
discharge  for  which  the  resonant  point  is  found  to  increase  in  magnitude 
and  move  to  a lower  frequency.  Increasing  the  slope  of  the  fan  map  tends 
to  flatten  the  pitch  response  and  also  shift  the  peak  response  to  a lower 
frequency . 

Modifications  in  fan  performance  may  be  expected  to  affect  heave 
response  more  than  pitch,  and  it  is  found  that  the  heave  response  to 
modifications  in  fan  discharge  is  roughly  opposite  that  of  the  pitch 
response.  As  the  flow  is  increased  the  magnitude  of  the  heave  response 
increases  at  frequencies  higher  than  the  resonant  frequency,  while  at 

i 

frequencies  below  the  resonance  range  there  is  almost  no  change  in 
response.  Finally,  there  appears  to  be  practically  no  effect  on  heave 
response  when  the  fan  slope  di/di  is  increased. 

SUBHARMONIC  OSCILLATIONS 

Oscillations  in  pitch  and  heave  at  frequencies  lower  than  the 
encounter  frequencies  are  predicted  by  the  analytical  model.  The 
presence  of  subharmonic  oscillations  is  quite  evident  in  the  pitch  and 
heave  time  histories  presented  in  Figure  17.  This  subharmonic  activity 
arises  from  the  nonlinear  character  of  the  analytical  model. 

Half  harmonic  oscillations  in  both  pitch  and  heave  begin  at 
nondimensional  or  reduced  frequencies  higher  than  7.0  for  a terrain  wave 
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amplitude  (a/L)  ratio  of  .02  and  a wave-to-cushlon-iength  ratio  of 
i.65.  The  magnitude  of  the  subharmonic  oscillations  increases  to  a 
maximum  at  a reduced  frequency  of  8.3  and  decreases  to  a negligible 
value  at  a reduced  frequency  of  12.  At  a smaller  wave  amplitude  ratio 
(a/L  • .01),  subharmonic  oscillations  do  not  appear  In  the  reduced 
frequency  range  considered  (0  to  15)  In  the  present  Investigation. 

However,  when  the  ambient  atmospheric  pressure  is  reduced  to  0.08  atm 
additional  subharmonics  occur,  ranging  from  the  1/3  to  1/5  subharmonics. 
These  additional  subharmonics  arc  attributed  to  compressibility  effects 
which  are  significant  at  reduced  pressures.  At  some  encounter  frequencies 
the  simulation  times  were  not  long  enough  to  reach  steady  periodic  motion 
in  pitch  or  heave.  This  occurred  at  the  points  In  the  encounter 
frequency  spectrum  where  the  response  oscillation  contained  more  than  one 
subharmonic  or  when  the  dominant  oscillations  were  undergoing  a transition 


from  one  predominant  subharmonic  to  another. 

CUSHION  COMPRESSIBILITY 

Compressibility  of  the  cushion  fluid  appears  to  have  a minor 
effect  on  model  scale  response  at  normal  ambient  (atmospheric)  pressure. 
However,  the  compressibility  effect  for  model  scale  predictions  at  a 
lower  ambient  pressure  (over-pressure),  corresponding  to  an  equivalent 
full-scale  operating  environment,  can  be  significant.  The  response  of 
the  craft  at  normal  atmospheric  pressure  is  compared  with  craft  response 
at  0.08  atm  In  Figures  l8  and  19.  At  the  lower  ambient  pressure  the 
effects  of  subharmonic  oscillations  are  significant  In  the  range  of 
wave-to-cushlon-length  ratios  between  1.5  and  2.5  (1.5  < A/L  < 2.5),  as 
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throughout  the  first  peak  for  both  heave  and  pitch  and  the  1/3  harmonic 
for  the  second  peak. 

Limited  analyticai  results  for  over-pressures  of  0.05  and  0.033 
atm  predict  a greatly  enhanced  response.  The  main  resonance  peak  magni- 
fications in  the  heave  and  pitch  transfer  functions  are  as  high  as  10  in 
certain  cases.  These  results  are  not  presented  here  since  the  goal  of 
the  examination  of  cushion  compressibility  Is  the  prediction  of  full- 
scale  ACV  response.  Other  components  of  the  ACV  dynamic  system  will  be 
modified  in  the  continuing  examination  of  scale  effects  in  order  to 
assess  the  exact  nature  of  the  effect  of  cushion  fluid  compressibility. 
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CONCLUSIONS 
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A successful  analytical  model  has  been  developed  to  predict  the 
nonlinear  pitch  and  heave  response  of  an  air-cushion-supported  craft. 
The  time-domain  analytical  model  is  appropriate  for  predicting  the 
dynamic  response  during  overland  operation.  Although  the  analytical 
model  is  quite  general,  the  motion  predictions  have  been  presented  for 
an  ACV  configuration  that  is  a suitable  representation  for  the  AALC 
Program  JEFF(B)  craft. 

The  predicted  results  obtained  for  both  pitch  and  heave  transfer 
functions  over  a complete  range  of  operating  conditions  are  compared 
with  available  experimental  data.  The  excellent  correlation  of  the 
results  establishes  the  accuracy  of  the  analytical  model  and  demon- 
strates its  usefulness  as  a research  and  development  tool.  To  this 
end,  various  parametric  sensitivity  studies  have  been  conducted  to  ex- 
plore the  effects  of  response  nonlinearity,  cushion  fluid  compressi- 
bility, and  variations  in  fan  performance.  The  result  of  the  linearity 
investigation  indicates  that  the  predicted  pitch  and  heave  transfer 
functions  have  regimes  of  linear  and  nonlinear  behavior  as  functions 
of  craft  speed,  amplitude,  and  wave  length  of  the  encountered  wave- 
field.  In  general,  the  nonlinearities  appear  to  dominate  for  small 
wave  slopes  (approaching  zero)  and  high  speeds. 

Cushion  fluid  compressibility  effects  appear  to  be  minor  for  a 
model-scale  craft.  However,  when  ambient  pressures  are  reduced  sig- 
nificantly in  order  to  effectively  model  a full-scale  cushion*  the 
effect  of  gas  compressibility  tends  to  enhance  the  pitch  and  heave 
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response.  Fairly  larqe  (approximately  15  per  cent)  parameter  vari- 
ations in  the  analytical  fan  model  can  be  tolerated  without  an  appre- 
ciable deqradation  in  predicted  craft  dynamic  performance.  This 
apparent  lack  of  sensitivity  to  fan  model  changes  indicates  that  fan 
system  parameter  optimization  to  decrease  power  requirements  may  be 
a possibi 1 i ty. 


RECOMMENDATIONS 


The  present  analytic  model  has  been  verified  on  a model  scale 
through  an  exceptionally  successful  correlation  between  numerical  and 
experimental  response  functions.  Based  on  these  results  it  is 
recommended  that  the  following  investigations  be  initiated. 

(1)  Extend  the  simulation  model  for  full-scale  JEFF(B)  char- 
acteristics and  perform  full-scale  predictions. 

(2)  Extend  the  model  to  include  the  cushion  pressure  gradient. 
Since  the  predicted  effects  of  gas  compressibility  proved 
more  significant  for  an  effective  full-scale  craft  than 
for  the  model,  it  is  expected  that  the  pressure  gradient 
associated  with  air  flow  in  the  cushion  may  have  a major 
influence  on  full-scale  craft  response. 

(3)  Extend  the  model  to  include  the  capability  of  predicting 
motion  in  all  six  degrees  of  freedom. 

(A)  Perform  additional  parameter  optimization  studies, 
especially  for  lift  system  performance. 
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^ (5)  Apply  the  predictive  analytical  model  to  the  AALC  Program 

JEFF(A)  craft. 

] • 

^ (6)  Examine  the  function  and  efficiency  of  proposed  vertical 

plane  control  devices  and  procedures.  Investigations  of 
this  type  are  meaningful  only  if  performed  on  a fully 
validated  mathematical  model. 
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Figure  2 - Bow  Skirt  Configuration  Employed  In  Analytica 
Mode  I 
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WAVE  LENGTH/CUSHION  LENGTH.  X/L 

Figure  3 - Correlation  of  Theoretical  Pitch  Transfer  Function 
with  Experimental  Data  as  a Function  of  Wave  Length 
for  Model  Speeds  of  k,  6,  and  8 Knots 
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WAVE  LENGTH/CUSHION  LENGTH,  X/L 

Figure  k - Correlation  of  Theoretical  Pitch  Transfer  Function 
with  Experimental  Data  as  a Function  of  Wave  Length 
for  Model  Speeds  of  10  and  12  Knots 
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WAVE  LENGTH/CUSHION  LENGTH,  X/L 


Figure  6 


- Correlation  of  Theoretical  Heave  Transfer  Function 
with  Experimental  Data  as  a Function  of  Wave  Length 
for  Model  Speeds  of  U,  6,  and  8 Knots 


HEAVE  TRANSFER  FUNCTION 
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Figure  8 


- Correlation  of  Theoretical  Heave  Transfer  Function 
with  Experimental  Data  as  a Function  of  Wave  Length 
for  Model  Speeds  of  13  and  15  Knots 


PITCH  TRANSFER  FUNCTION 


Figure  9 ~ Response  Linearity:  Pitch  Transfer  Function  as  a 


Function  of  Wave  Slope  for  Model  Speeds  of  11 , 13.5, 
16,  and  18  Knots 
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Figure  10  - Response  Linearity:  Heave 

Function  of  Wave  Slope  for 
16,  and  18  Knots 
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NONDIMENSIONAL  ENCOUNTER  FREQUENCY  c^s/Ug 

Figure  11  - Response  Linearity:  Pitch  Transfer  Function  as  a 

Function  of  Nondimensional  Encounter  Frequency  for 
X/L  * 1.65 
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Figure  12  - Response  Linearity:  Heave  Transfer  Function  as  a 

Function  of  Nond i mens  Iona  1 Encounter  Frequency  for 
X/L  = 1.65 
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FLOW  COEFFICIENT  RATIO  <I>AI> 


Figure  14  - Fan  Performance  Curves  Used  In  Parametric  Study  of 
Lift  Fan  Characteristics 
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Figure  16  - Parametric  Study  of  Lift  Fan  Characteristics: 

Heave  Transfer  Function  as  a Function  of  Encounter 
Frequency  for  Various  Fan  Maps  with  X/L  = 1.65 
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Time-Dependent  Pitch  and  Heave  Simulation  Responses 
for  Sinusoidal  Terrain  at  a Model  Speed  of  18  Knots 


HEAVE  TRANSFER  FUNCTION 


Transfer  Function  as  a Function  of  Wave  Length  for  a • 

Model  Speed  of  15  Knots 
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